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(57) Abstract 



A transmitter (10) provides an output signal indicative of mass flow rate of fluid through a conduit (12). The tran.smitter (10) 
includes a temperature sensor (34) providing a tcmpcramre signal indicative of fluid temperature. A static pressure sensor (30) provides a 
static pressure signal indicative of static pressure in the conduit. A differential producer (32) provides a differential pressure signal The 
transmitter (10) also includes a controller (66) which provides the output signal indicative of mass flow of the fluid through the conduit. 
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TRANSMITTER FOR PROVIDING A SIGNAL INDICATIVE OF FLOW THROUGH A 
DIFFERENTIAL TRANSDUCER USING A SIMPLIFIED PROCESS 



BACKGROUND OF THE INVENTION 
5 The presenc invention deals with a transiritter 

in the process control industry. More particularly, the 
present invention deals with a simplified process, used 
in a transmitter, for providing an output signal 
indicative of flow through a differential producer. 

10 Transmitters which sense various 

characteristics of fluid flowing through a conduit are 
known. Such transmitters typically sense and measure 
differential pressure, line pressure (or static 
pressure) and temperature of the process fluid. Such 

15 transmitters are typically mounted in the field of a 
refinery, or other process control industry 
installation. The field mounted transmitters are 
subject to significant constraints on power consumption. 
Such transmitters commonly provide an output in the form 

20 of a current representative of the variable being 
sensed. The magnitude of the current varies between 
4-20 mA as a function of the sensed process variable. 
Therefore, the current available to operate the 
transmitter is less than 4 mA. 

25 One way in which flow computation is done in 

industries such as the process control industry and the 
petroleum industry is through the use of dedicated flow 
computers. Such devices either use separate pressure, 
differential pressure and temperature transmitters or 

30 have sensing mechanisms housed in large enclosures. 
These devices are generally large and consume more power 
than 4 mA. Addi ticr.ally , they are often limited to use 
in specialized applications such as the monicorinq cf 
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hydrocarbons for custody transfer or at wellheads to 
monitor the output of gas or oil wells. 

Another way in which flow computation is done 
is through the use of local control systems, often 
5 called programmable loop controllers (PLC). PLC's 
typically receive inputs from separate pressure, 
differential pressure and temperature transmitters and 
compute the flow based on these inputs. Such devices 
are often performing additional local control tasks such 
.0 as the calculation of other variables required in the 
control of the plant or the monitoring of process 
variables for alarm purposes. The calculation of flow 
in these devices requires programming by the user. 

A third way in which flow computation is done 
.5 is through the use of large computers which control 
entire plants, often called distributed control systems 
(DCS). DCS's typically perform a wide range of tasks 
ranging from receiving inputs from field-based 
transmitters to computing the intermediate process 
0 variables such as flow or level, to sending positioning 
signals to final control elements such as valves, to 
performing the monitoring and alarm functions within the 
plant . Because of the wide range of tasks required and 
the typically high cost of DCS input/output capability, 
5 memory and computational time, it is common to do a flow 
computation that is not compensated for all of the 
effects due to changing process conditions. 

One common means of measuring flow rate in the 
process control industry is to measure the pressure drop 
0 across a fixed restriction in the pipe, often referred 
to as a differential producer or primary element. The 
general equation for calculating flow rate through a 
differential producer can be written as: 



wo 97/04288 PCTAJS96/1 1 51 5 



-3- 



Equation 1 



0- NC^EY^d^ yfph 



where 



Q = Mass flow rate (mass/unit time) 
5 N = Units conversion factor (units vary) 

Cd = Discharge coefficient (dimensionless) 
E = Velocity of approach factor 

(dimensionless) 
Yi = Gas expansion factor (dimensionless) 
10 d = Bore of differential producer (length) 

p = Fluid density (mass/unit volume) 
h = Differential pressure (force/unit area) 
Of the terms in this expression, only the 
units conversion factor, which is a constant, is simple 
15 to calculate. The other terms are expressed by 
equations that range from relatively simple to very 
complex. Some of the expressions contain many terms and 
require the raising of numbers to non-integer powers. 
This is a computationally intensive operation. 
20 In addition, it is desirable to have the 

transmitter operate compatibly with as many types of 
differential producers as possible. Implementing all of 
the calculations and equations needed for the 
conventional flow equation in order to determine flow 
25 based on the output of one differential producer (much 
less a plurality of different types of differential 
producers) requires computations which can only be 
reasonably performed by a processor which has a high 
calculation speed and which is quite powerful, 
30 Operation of such a processor results in increased power 
consumption and memory requirements in the transmitter. 
This is highly undesirable given the 4 mA power 
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constraint or conventional transmitters. Therefore, 
current transmitter-based microprocessors, given the 
above power and memory constraints, simply do not have 
the capability of performing the calculations in any 
5 reasonable time period. 

There has been some work done in obtaining a 
simplified discharge coefficient equation. However, 
this is only one small part of the flow equation. Even 
assuming the discharge coefficient is extremely 
10 simplified, implementing the flow equation accurately is 
still very difficult given the constraints on current 
transmitter-based microprocessors . 

Other attempts have been made to simplify the 
entire flow equation. However, in order to make the 
15 flow equation simple enough that it can be implemented 
in transmitter-based microprocessors, the simplified 
flow ecjuations are simply not very accurate. For 
example, some such simplified flow equations do not 
account for the discharge coefficient. Others do not 
20 account for compressibility, or viscosity effects. 

Therefore , common transmitter -based 
microprocessors which are powered by the 4-20 mA loop 
simply do not accurately calculate flow. Rather, they 
provide outputs indicative of differential pressure 
25 across the orifice plate, static line pressure, and 
temperature. These variables are provided to a flow 
computer in a control room as mentioned above, which, in 
turn, calculates flow. This is a significant processing 
burden on the flow computer. 
30 SUMMARY OF THE INVENTION 

A transmitter provides an output signal 
indicative of mass flow rate of fluid through a conduit. 
The transmitter includes a temperature sensor providing 
a temperature signal indicative of fluid temperature. 
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A static pressure sensor provides a static pressure 
signal indicative of static pressure in the conduit. A 
differential producer provides a differential pressure 
signal. The transmitter also includes a controlleir 
5 which provides the output signal indicative of mass flow 
of the fluid through the conduit based on a plurality of 
simplified equations . 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 shows a transmitter according to the 
10 present invention connected to a pipe which conducts 
fluid therethrough. 

FIG. 2 is a block diagram, in partial 
schematic form of the transmitter according to the 
present invention. 
15 FIGS. 3A-3C graphically illustrate curve fit 

accuracy for the discharge coefficient used by the 
system according to the present invention. 

FIGS. 4A and 4B graphically illustrate curve 
fit accuracy of viscosity used according to the present 
20 invention. 

FIG. 5 illustrates the curve fit accuracy of 
the term Ed^ used according to the present invention. 

FIG. 6 graphically illustrates the curve fit 
accuracy of the gas expansion factor used according to 
25 the present invention. 

FIGS. 7A and 7B graphically illustrate the 
curve fit accuracy of fluid density for liquid used 
according to the present invention. 

FIGS. 8A and SB graphically illustrate curve 
30 fit accuracy of fluid density for gas used according to 
the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 



FIG. 1 is an illustration of a transmitter 10 
according to the present invention. Transmitter 10 is 
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coupled to a pipe 12 through pipe fitting or flange 14. 
Pipe 12 conducts flow of a fluid, either a gas or a 
liquid, in the direction indicated by arrow 16. 

Transmitter 10 includes transmitter 
5 electronics module 18 and sensor module 22 which 
collectively house a transmitter more fully illustrated 
in FIG. 2. Transmitter electronics module 18 also 
preferably includes a boss 20 for accepting an input 
from a resistive temperature device (RTD) , preferably a 

10 100 ohm RTD which is typically inserted directly into 
the pipe or into a thermowell which is inserted into the 
pipe to measure the fluid temperature. The wires from 
the RTD are connected to one side of a terminal block in 
a temperature sensor housing 24, To the other side of 

15 the terminal block are connected wires which run through 
an electricl conduit 26 and are coupled to boss 20. 

Sensor module 22 includes a differential 
pressure sensor and an absolute pressure sensor. The 
differential pressure sensor and absolute pressure 

20 sensor provide pressure signals to conditioning and 
digitizing circuitry, and to a linearizing and 
compensating circuit. The compensated, linearized and 
digitized signals are provided to the electronics module 
18. The electronics module 18 in transmitter 10 

25 provides an output signal indicative of process 
conditions of the process fluid flowing through pipe 12 
to a remote location, by a 4-20 mA two- wire loop 
preferably formed using twisted pair conductors, through 
flexible conduit 28. In the preferred embodiment, 

30 transmitter 10 provides signals which are indicative of 
the three process variables (temperature, static 
pressure, and differential pressure) according to the 
HAJIT® or Fieldbus Standards. Further, in accordance 
with the present invention, transmitter 10 also provides 
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an output signal indicative of flow. The method of 
determining flow according to the present invention is 
significantly simplified over prior methods allowing the 
microprocessor in the electironics tnodule of transmitter 
10 to calculate flow without exceeding the power 
constraints on the microprocessor, and at acceptably 
fast update times. 

FIG. 2 is a more detailed block diagram of 
sensor module 22 and electronics module 18 of 
transmitter 10. Sensor module 22 includes a strain 
gauge pressure sensor 30, differential pressure sensor 
32 and temperature sensor 34. Strain gauge sensor 30 
senses the line pressure (or static pressure) of fluid 
flowing through conduit 12. Differential pressure 
sensor 32 is preferably formed as a metal cell 
capacitance-based differential pressure sensor which 
senses the differential pressure across an orifice in 
conduit 12. Temperature sensor 34, as discussed above, 
is preferably a 100 ohm RTD sensor which senses a 
process temperature of fluid in pipe 12. While, in FIG. 
1, sensor 34 and sensor housing 24 are shown downstream 
of transmitter 10, this is but one preferred embodiment, 
and any suitable placement of temperature sensor 34 is 
contemplated. 

Sensor module 22 also preferably includes an 
analog electronics portion 36, and a sensor processor 
electronics portion 38. Electronics module 18 includes 
output electronics portion 40. Analog electronics 
portion 36 in sensor module 22 includes signal 
conditioning and power supply filtering circuitry 42, 
analog-to-digital (A/D) circuitry 44, and PRT 46. The 
analog signals received from sensors 30, 32 and 34 are 
provided to analog signal conditioning and power supply 
filtering circuitry 42. The analog signals are 
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conditioned (such as amplified) and the conditioned 
analog signals are provided to A/D converter circuitry 
44. 

In a preferred embodiment, A/D converter 
5 circuitry 44 includes a plurality of voltage- to-digital 
converters, or capacitance- to-digital converters, or 
both (as appropriate) which digitize the analog input 
signals. Such converters are preferably constructed 
according to the teachings of U.S. Patent Nos . 

10 4,878,012; 5,083,091; 5,119,033 and 5 , 155 , 455 ; assigned 
to the same assignee as the present invention. In the 
embodiment shown in FIG. 2, three voltage- to-digital 
converters 48, 50 and 52, and one capacitance- to-digital 
converter 54 are shown'. The voltage -to-digital 

15 converters 48 and 50 are used to convert the signals 
from sensors 30 and 34 into digital signals. The 
capacitance -to-digital converter 54 is used to convert 
the signal from capacitive pressure sensor 32 to a 
digital signal. 

20 PRT 4 6 is preferably formed as a low cost, 

silicon-based PRT positioned proximate pressure sensors 
30 and 32. PRT 46 provides a temperature signal 
indicative of the temperature proximate sensors 30 and 
32. This temperature signal is provided to voltage-to- 

25 digital converter 52 where it is digitized. This 
digitized signal is then used to compensate the 
differential and line pressure signals for temperature 
variations. Analog signal conditioning and power supply 
filtering circuitry 42, the A/D converters 44 and PRT 46 

30 are all preferably physically located proximate to, or 
on, a single circuit board housed in transmitter 10. 

Once the analog signals are digitized by A/D 
converters 44, the digitized signals are provided to 
sensor processor electronics portion 38 as four 
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respective sixteen bit wide outputs on any suitable 
connection or bus 56. 

Sensor processor electronics portion 38 
preferably includes a microprocessor 58, clock circuitry 
5 60 and memory (preferably electrically erasable 
programmable read only memory, EEPROM) 62. 
Microprocessor 58 compensates and linearizes the process 
variables received from analog electronics portion 36 
for various sources of errors and non-linearity. For 

10 instance, during manufacture of transmitter 10, pressure 
sensors 30 and 32 are individually characterized over 
temperature and pressure ranges, and appropriate 
correction constants are determined. These correction 
constants are stored in EEPROM 62. During operation of 

15 transmitter 10, the constants in EEPROM 62 are retrieved 
by microprocessor 58 and are used by microprocessor 58 
in calculating polynomials which are, in turn, used to 
compensate the digitized differential pressure and 
static pressure signals. 

20 Clock circuitry 60 is provided in sensor 

processor electronics portion 38 and provides clock 
signals to microprocessor 58, A/D circuits 44 and to 
other electronics as appropriate, in order to accomplish 
the desired operations. It should also be noted that 

25 the functionality of portions 36 and 38 can be combined 
into a single integrated circuit chip through 
application specific integrated circuit (ASIC) 
technology. 

After the analog signals from sensors 30, 32 
30 and 34 are digitized, compensated and corrected, the 
process variable signals are provided over a serial 
peripheral interface (SPI) bus 64 to output electronics 
portion 40 in electronics module 18. SPI bus 64 
preferably includes power signals, two hand shaking 



wo 97/04288 



PCTAJS96/11515 



-10- 

signals and the three signals necessary for typical SPI 
signaling. 

Output electronics module 4 0 preferably 
incliades Tnicroprocessor 66, non-volatile memory 6B, 

5 voltage regulator 72, modulator circuit 74, HART® 
protocol receiver 76 and loop current controller 78. In 
addition, output electronics portion 40 may optionally 
be coupled to a battery back-up circuit which provides 
battery power to the output electronics in case of 

0 failure of the power provided over the two-wire loop. 

Microprocessor 66 receives the digitized, 
compensated process variables over SPI bus 64 . In 
response, and as will be described in greater detail 
later in this specification, microprocessor 66 

5 calculates the mass flow of fluid flowing through pipe 
12 based on the process variables received over bus 64 . 
This information is stored in non- volatile memory 6 8 
which, preferably, is suitable for storing up to 35 days 
worth of mass flow data, 

0 When requested, microprocessor 66 configures 

output electronics 4 0 to provide the mass flow data 
stored in non-volatile memory 68 over two-wire loop 82. 
Therefore, output electronics 40 is coupled at positive 
and negative terminals 84 and 86 to loop 82 which 

5 includes controller 88 (modeled as a power supply and a 
resistor). In the preferred embodiment, output 

electronics 40 communicates over two-wire loop 82 
according to a HART® communications protocol, wherein 
controller 88 is configured as a master and transmitter 

0 10 is configured as a slave. Other communications 
protocols common to the process control industry may be 
used, with appropriate modifications to the code used 
with microprocessor 66 and to the encoding circuitry. 
Communication using the HART® protocol is accomplished 
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by utilizing HART® receiver 76. HART® receiver 76 
extracts digital signals received over loop 82 from 
controller 88 and provides the digital signals to 
circuit 74 which, in turn, demodulates the signals 
5 according to the HART® protocol and provides them to 
microprocessor 66 . 

Circuit 74 receives digital signals (which are 
to be sent over loop 82) from microprocessor 66. 
Circuit 74 converts the digital signals into analog 

10 signals, modulates them for transmission, and provides 
the modulated signals to circuit 76. Circuit 74 
preferably includes a Bell 22 compatible modem. The 
loop current control circuit 78 receives an analog 
voltage signal from a D/A converter in circuit 74. In 

15 response, loop current control circuit 78 provides a 4- 
20 mA output representative of the particular 
information being transmitted by microprocessor 66 over 
loop 82 (such as one of the process variables, or the 
calculated flow) . 

20 Also, voltage regulator 72 preferably provides 

3.5 volts and other suitable reference voltages to 
output electronics circuity 40, sensor processor 
electronics 38, and analog electronics 36. 

In order to calculate flow through a 

25 differential producer (such as an orifice plate) 
information is required about three things. Information 
is required about the process conditions, about the 
geometry of the differential producer and about the 
physical properties of the fluid. Information about the 

30 process conditions is obtained from sensor signals, such 
as the signals from sensors 30, 32 and 34. Information 
regarding the geometry of the differential producer and 
the physical properties of the fluid are provided by the 
user. 
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Flow through a differential producer is 
conventionally calculated by utilizing the equation set 
out as Equation 1 above. Flow is typically calculated 
in mass units, but can be expressed in volumetric units 
if required. The choice of units determines the value 
of the units conversion factor, N. 

The discharge coefficient, C^, is a 
dimensionless, empirical factor which corrects 
theoretical flow for the influence of the velocity 
profile of the fluid in the pipe, the assumption of zero 
energy loss in the pipe, and the location of pressure 
taps. Cd is related to the geometry of the differential 
producer and can be expressed as a seemingly simple 
relationship in the following form: 

Equation 2 

where the Reynolds number 

Re„" = ^ ; 

\iD 

C„ = the discharge coefficient at infinite 
20 Reynolds number; 

b = a known Reynolds number correction term; 
n = a known exponent term; and 
/i = the fluid viscosity. 

This relationship varies for different types 
25 of differential producers, the location of the pressure 
taps on such producers, and the beta ratio. Typical 
equations defining C^, and the other above terms have a 
wide range of complexity and are set out in Table 1. 
The calculation for associated with an orifice plate- 



10 
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type differential producer is the most common in the 
industry. 

The velocity of approach factor, is a 

geometrical term and relates the fluid velocity in the 
throat of the differential producer to that in the 
remainder of the pipe. The velocity of approach factor 
is a function of temperature as follows: 

Equation 3 




where, for an orifice meter, 

Equation 4 

d,[l + a, (T- T,) ] 



P- 



dj. = orifice diameter at reference temperature 

T • 

Dj. = meter tube diameter at reference 
temperature T^; 

= thermal expansion coefficient of the 
orifice plate; and 

a2 = thermal expansion coefficient of a meter 

tube . 

The gas expansion factor Yj is a dimensionless 
factor which is related to geometry, the physical 
properties of the fluid and the process conditions. The 
gas expansion factor accounts for density changes as the 
fluid passes through a differential producer. The gas 
expansion factor for primary elements with abrupt 
changes in diameter, such as orifice meters, is given by 
the following empirical relationship: 
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° 1 - (.41 + .35p^) 



27 .73PK 



Equation 5 



where h = differential pressure in inches of 
water at 68^F; 

P = upstream pressure in psia; and 

K = isentropic exponent of the gas. 

The adiabatic gas expansion factor for 
contoured elements is described as follows: 

Equation 6 



^1 = 



10 



where 



(1 - P^) [K/K-D] {P.,/P^)^/^[l - (P^/P,) 



[1 - pMPj/Pi) (I-P2/P,) 



1/2 



Equation 7 



^ - 1 



27 ,73P 



K = isentropic exponent of the gas. 
15 The value of is 1.0 for liquids. 

The bore of the differential producer, d, is 
related to geometry and is a function of temperature as 
follows : 

Equation 8 

d-d^[l + a,(r- rj] 

20 

The differential pressure factor, h, is 
measured by a conventional differential pressure sensor. 

The fluid density factor p is described in 
mass per unit volume and is a physical property of the 
25 fluid. For typical process control applications, the 
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density of liquids is a function of temperature only. 
It can be described by expressions such as the PTB 
equation for the density of water: 

Equation S 

p « A BT + CT^ + DT^ + ET"^ + FT^ 

where A-F are constants, or a generic 
expression given by the American Institute of Chemical 
Engineers (AIChE) : 

Equation 10 

aM 



-T/C) ^ 



10 

Where a-d are fluid dependent constants and M 
is the molecular weight. 

Gas density is a function of absolute pressure 
and absolute temperature given by the real gas law: 
15 Equation 11 

P 



P " 



nZRJ 



where Z = the compressibility factor; 

Rq = universal gas constant ; and 

n = number of moles . 
20 Gas density and compressibility factors are 

calculated using equations of state. Some equations of 
state, such as AGA8, the ASME steam equation and MBWR, 
are useful for single fluids or a restricted number of 
fluids. Others, such as Redlich-Kwong or AIChE 
25 equations of state are generic in nature and can be used 
for a large number of fluids. The AIChE equation is as 
follows : 
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Equation 12 



_ 1 _ 


1 


1 


4P 


1/2 


~ 2S~ 


2 




BRT 





where 

Equation 13 




5 

where a-e are fluid dependent constants; and 
M = the molecular weight of the fluid. 
Implementing the flow calculation using 
equations 1-13 set out above would yield a highly 
10 accurate result. However, the constraints of power 
consumption, calculation speed and memory requirements 
make the implementation of the full equations beyond the 
capability of currently available transmitter based 
microprocessors. Therefore, the transmitter of the 
15 present invention calculates flow based on a number of 
simplified equations, while retaining a high degree of 
accuracy in the flow calculation. 

The dependencies related to the discharge 
coefficient are as follows: 
20 ()S, Rer>) ; 

Reo (Q, m) where /z is the viscosity 

of the fluid; and 
^l (T) 

Using the AIChE equation for liquids: 
25 Equation 14 

M = exp (a + b/T + c Ln(T) + d T^) ; 
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and, using AIChE equation for gases: 

Equation 15 



10 



15 



1 + c/T + d/T^ 



According to the present invention, the 
discharge coefficient equation is simplified by 

approximating /i"* by a polynomial in T or 1/T. 
Preferably, this approximation is done using a third 
degree polynomial equation. Also, is approximated 
using a sixth degree polynomial equation in 

1 



or 



LniReJ 



It has been observed that better accuracy is 
obtained using the polynomial for Q with the term 

1 

LniRe^) 



being the independent variable, but this also increases 
the calculation time. Therefore, this can be used, or 
the other polynomial can be used, depending upon the 
degree of accuracy desired. 
20 FIGS. 3A, 3B and 3C are examples of curve fit 

accuracy of the discharge coefficient using the above 
equations. FIG. 3A is a graph of the discharge 
coefficient curve fit error versus the Reynolds number 
for an ASME flange tap orifice meter having a diameter 
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in excess of 2.3 inches. This graph was obtained by 
doing a sixth degree fit in 

1 



as follows; 



\fRei 
1 

/ReA ' '/Re, 



y/Re, 



Equation 16 



y/Re 



DJJJ 



and using an approximation of viscosity as follows: 

Equation 17 



10 FIG. 3B graphically illustrates the discharge 

coefficient curve fit error plotted against Reynolds 
number for an ASME corner tap orifice meter using a 
sixth degree fit in 

1 



15 FIG. 3C graphically illustrates the discharge 

coefficient curve fit error against Reynolds number for 
an ASME long radius nozzle using the sixth degree fit in 



1 
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FIGS. 3A-3C illustrate that the curve fit 
approach approximates the discharge coefficient to 
better than +/-0.005%. Similar results are obtained for 
other differential producers. 
5 FIGS. 4A and 4B are examples of curve fit 

accuracy obtained for viscosity. FIG. 4A graphically 
illustrates curve fit accuracy for viscosity versus 
temperature using a third degree polynomial fit in l/T. 
FIG. 4B illustrates the curve fit accuracy for viscosity 

10 versus temperature using a third degree polynomial fit 
in l/T. FIG. 4A is based on water and FIG. 4B is 
calculated for air. It is seen that the curve fit 
approach approximates the viscosity of air to better 
than +/- 0.001% and the viscosity of water to better 

15 than +/- 0.2%. A polynomial fit of a higher degree in 
l/T, such as 4 or 5, would improve the accuracy of the 
fit for water. Because the discharge coefficient, C^, 
is weakly dependent on Reynolds number and, thus, 
viscosity, the accuracy provided using a third degree 

20 polynomial fit in l/T is acceptable and the added 
computational complexity of a higher degree polynomial 
approximation is not necessary. Similar results are 
obtained for other liquids and gases. 

The dependencies related to the velocity of 

25 approach factor, E, and the bore of the differential 
producer, d, are as follows: 

E(T) , and dMT) . 

30 The method of the present invention simplifies 

the Ed^ calculation by grouping E and d^ together and 
approximating the product of Ed^ by a polynomial in T or 
l/T. This polynomial is preferably a second degree 
polynomial . 
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FIG. 5 is an example of the curve fit accuracy 
of the Ed^ term. FIG. 5 graphically illustrates the 
accuracy of this term plotted against temperature using 
a second degree polynomial in T as follows: 

Equation 18 



£■^2 ^ Cn - c 



FIG. 5 illustrates that the curve fit approach 
approximates the Ed^ term to better than +/-b. 00002%. 

The dependencies of the gas expansion factor, 
10 Yi, are as follows: 

r,(p,K, |) ;ancl 



Simplifying the gas expansion factor 
calculation is accomplished by ignoring the dependency 

15 on T. The term is approximated using a polynomial 
equation in h/P where h is the differential pressure and 
P is the static pressure. Preferably, this polynomial 
is a second degree polynomial. For an orifice, a linear 
relationship exists between Yj and h/P. 

20 FIG. 6 is an example of curve fit accuracy of 

Yj versus temperature using a second degree polynomial 
fit in h/P as follows: 

Equation 19 

25 The curve is illustrated for a contoured 

element differential producer. FIG. 6 illustrates that 
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the term for contoured elements is accurate, using 
the system according to the present invention, to better 
than +/-0.002% for all beta ratios. Accuracy is better 
than +/-O.0OO5% for beta ratios less than 0.6. Similar 
5 results are obtained for the square edged orifice. 

Dependencies related to the fluid density for 
liquid and gas are as follows: 

Pug(T) ; and 

Pa.siP''^)' specifically p^,, (-^]^ 

The fluid density calculation for liquid is 
simplified according to the present invention by 
providing two levels of curve fit. The term /pJJ^ is 
approximated by a polynomial in T or 1/T. Preferably, 
this is a third degree polynomial and is provided as a 
default equation for a lower accuracy fit as follows: 

Equation 20 

The same term is also preferably approximated by a 
polynomial in l/T using a fifth degree polynomial as a 
20 higher accuracy fit for broader operating ranges of 
temperature . 

Simplifying the calculation for fluid density 
for gas is accomplished by, again providing two levels 
of curve fit. Fitting a curve to — and not p^^^ 

25 improves the curve fit accuracy, reduces calculation 
time, and improves the simplified flow equation 

accuracy. According to the present invention, the term — 

/z 

is approximated by a polynomial in P and 1/T. In the 



10 



15 
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preferred embodiment, the default polynomial is a 3X2 
polynomial and is used for a lower accuracy fit. 
However, the term — can also be approximated by a 

polynomial in P and 1/T using an 8X6 polynomial for 
higher accuracy fits, and for broader operating ranges 
of both P and T. The preferred simplified equation for 
fluid density for all gases is as follows: 

Equation 21 



v/p 



144A/„" 


.5 


P' 


R 




T 



10 FIG. 7A graphically illustrates an example of 



curve fit accuracy for yfpJJ^ for water versus 
temperature using the third degree polynomial fit in 
l/T. FIG. 7B graphically illustrates curve fit accuracy 
density of acrylonitrile versus temperature. In both 
15 cases, the temperature range is 50°F to llO^F. FIGS- 7A 
and 7B illustrate that the curve fit approach 



approximates v/Pjiq better than +/- 0.0002% for these 
two liquids and the selected temperature range. Similar 
results are obtained for other liquids and other 
20 temperature ranges, 

FIGS. 8A and SB illustrate examples of curve 

fit accuracy for — for two fluids and pressure 

/i 

temperature ranges. FIG. 8A illustrates the curve fit 
accuracy using the 3X2 polynomial fit for carbon dioxide 
25 gas. The pressure and temperature ranges are 15 psia to 
115 psia and 60°F to 140°F. The results show that the 

curve fit approach accurately approximates — to 

v/z 

better than +/- 0.0015%. FIG. 8B illustrates the curve 
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fit accuracy using the 3X2 polynomial fit for ethylene 
gas. The pressure and temperature ranges are 75 psia to 
265 psia and GO^F to 140®F. The results show that the 

curve fit approach accurately approximates — to 

5 better then +/- 0.005%. As these results indicate, the 
accuracy of the curve fit approximation varies, as the 
fluid is changed and as the operating ranges of pressure 
and/or temperature change. When the operating ranges of 
pressure and/or temperature result in unacceptable 

10 approximations by using a 3X2 polynomial, an 8X6 
polynomial will improve the results to levels similar to 
those indicated in FIGS. 8A and 8B. 

In sum, the classic flow calculation given by 
Equation 1 above, is simplified according to the present 

15 invention as follows: 

Equation 22 

For gases this equation can be rewritten as: 

1 



0 - KN [CJ lEd'] [yj 
where 



hp 
T 



144^ 



20 M = molecular weight of the gas; 

R = gas constant; and 

P, h, T are in units of psia, inches of water 
and degrees Rankine, respectively. For liquids, the 
equation can be rewritten as: 

Q - N ICJ [Ed^] [yj l/p]y/h 
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where the bracketed terms are curve fit approximations. 
By simplifying the flow equation as set out above, the 
transmitter based microprocessor 66 is capable of 
updating the flow calculation each time it receives 
5 updated sensor information by bus 64. In the event that 
one or more of the curve fit approximations have not 
been completely calculated the previous value is used in 
the flow calculation. 

The effect of variations in the process 

10 variables has a direct affect on the flow calculation by 
virtue of their appearance in the flow equation. They 
have a smaller effect on the curve fit cerms. Thus, by 
using the newly updated process variable information and 
the most recently calculated curve fit approximations, 

15 the result is a flow calculation that is both fast and 
accurate. Having newly calculated flow terms at such an 
expedient update rate allows transmitter 10 to exploit 
fast digital communication protocols. 

Also, by simplifying the flow calculation as 

20 set out above, microprocessor 66 performs the same 
calculations regardless of the type of differential 
producer used, regardless of the beta ratio used, and 
regardless of whether the user requires a simplified or 
fully compensated flow. 

25 It should also be noted that the curve fit 

coefficients are easily calculable by the user using 
known techniques. These coefficients are simply stored 
in memory associated with microprocessor 66 and used in 
performing the desired calculations. 

30 These simplifications allow transmitter 10 to 

actually calculate flow in a highly accurate manner. 
Rather than requiring the transmitter to simply transmit 
the process variables back to a control room, and have 
a flow computer in the control room or installation 
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calculate flow, the transmitter according to the present 
invention is capable of not only providing the process 
variables, but also providing a flow calculation to the 
control room. This relieves the processing overhead on 
5 the flow computer or other processor in the control 
room, yet does not over burden the transmitter-based 
microprocessor, or require the transmitter-based 
microprocessor to use energy which exceeds that 
available to it. 
10 Although the present invention has been 

described with reference to preferred embodiments, 
workers skilled in the art will recognize that changes 
may be made in form and detail without departing from 
the spirit and scope of the invention. 
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WHAT IS CLAIMED IS: 

1. A transmitter for providing an output signal 
indicative of mass flow rate of fluid through a conduit, 
the transmitter comprising: 

a temperature sensor providing a temperature 
signal indicative of fluid temperature; 

a static pressure sensor providing a static 
pressure signal indicative of static 
pressure in the conduit; 

a differential producer providing a 
differential pressure signal; and 

a microcomputing circuit, coupled to the 
temperature sensor, the static pressure 
sensor, and the differential producer, 
to receive the temperature signal, the 
static pressure signal and the 
differential pressure signal, and 
providing an output signal indicative of 
flow of the fluid through the conduit 
based on a plurality of polynomial 
equations using polynomial curve fits 
with at least one of temperature, static 
pressure and differential pressure being 
an independent variable in each of the 
polynomial equations . 

2. The transmitter of claim 1 wherein the 
microcomputing circuit includes: 

a first microprocessor coupled to the 
temperature sensor, static pressure 
sensor and differential producer, and 
corrects the static pressure signal, 
differential producer signal and 
temperature signal for non-linearities 
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and provides corrected output signals; 
and 

a second microprocessor, coupled to the first 
microprocessor, for calculating flow 
based on the corrected output signals. 
3. A 4-20 mA process control transmitter coupled 

to a conduit conducting a fluid therethrough, the 
transmitter comprising: 

a first pressure sensor sensing line pressure 
in the conduit and providing a line 
pressure signal indicative of the line 
pressure; 

a second pressure sensor sensing differential 
pressure across an orifice in the 
conduit and providing a differential 
pressure signal indicative of the 
differential pressure; 
a temperature sensor sensing temperature of 
the fluid and providing a temperature 
signal indicative of the temperature of 
the fluid; and 
a microcomputing circuit, coupled to the 
first and second pressure sensors and 
the temperature sensor and powered by a 
4-20 mA loop, calculating flow of the 
fluid through the conduit based on the 
line pressure signal, the differential 
pressure signal and the temperature 
signal and providing an output signal 
indicative of the flow. 
4. The transmitter of claim 3 wherein the 

controller is configured to calculate flow based on at 
least one polynomial equation using a polynomial curve 
fit with at lease one of temperature, static pressure 
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and differential pressure being an independent variable 
in the polynomial equation. 

5. The transmitter of claim 4 wherein the 
controller calculates flow, Q, according to an equation 
generally of the form: 

0-N[C^] [Ed^] [yj [x/p]v/7i 

wherein N is a units conversion factor, is a 

discharge coefficient, E is a velocity of approach 
factor, d is a bore diameter of the second pressure 
sensor, is a gas expansion factor, p is fluid density 
at operating temperature and pressure, and h is the 
differential pressure, and wherein at least one of C^, 
E, d^, y^, and p are approximated by the polynomial 
equation curve fit. 

6. The transmitter of claim 5 wherein the 
controller calculates flow based on a plurality of 
polynomial equations using polynomial curve fits to 
approximate a plurality of C^, E, d^, Y^, and p with at 
least one of the temperature, the static pressure and 
the differential pressure being an independent variable 
in the polynomial equations. 

7. The transmitter of claim 5 wherein the 
controller calculates flow based on an approximation of 
Cd according to a polynomial equation having the form: 
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8 . The transmitter of claim 7 wherein is 

calculated as : 



yfRe 



A * — ( 



dJ/I) 



9. The transmitter of claim 5 wherein the 

controller calculates flow based on an approximation of 
Ed^ according to a polynomial equation having the form: 



10. The transmitter of claim 9 where Ed'' is 

calculated as : 



11. The transmitter of claim 5 wherein the 

controller calculates flow based on an approximation of 
Yi according to a polynomial equation having the form: 



12. The transmitter of claim 11 wherein is 

calculated as: 



13. The transmitter of claim 5 wherein the 

controller calculates flow based on an approximation of 
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p for liquid according to a polynomial equation having 
the form: 



14. The transmitter of claim 13 wherein p for 

liquid is calculated as: 



15. The transmitter of claim 5 wherein p for a gas 

is calculated as: 



v/p - 



14 4W^,' 


.5 


■ P' 


L R \ 




T 



'LLP' (iff. 

i-0 j-0 \ ^ / 



16 . The transmitter of claim 15 wherein p for gas 

is calculated as: 





.5 


■ P 


R 







.5 



17. The transmitter of claim 5 wherein is 

calculated using an equation generally in the form: 
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18. The transmitter of claim 16 wherein is 
calculated as: 

19. The transmitter of claim 5 wherein the term 
Ed^Yi is calculated using an equation in the form: 




20. A method of providing an indication of flow of 

fluid through a conduit using a process control 
transmitter powered by a 4 -20mA control loop, 
comprising: 

sensing static pressure and providing a 
static pressure signal indicative of the 
static pressure; 
sensing differential pressure and providing a 
differential pressure signal indicative 
of the differential pressures- 
sensing temperature and providing a 
temperature signal indicative of the 
temperature; and 
calculating flow of the fluid through the 
conduit based on the static pressure 
signal, the differential pressure signal 
and the temperature signal and providing 
an output signal indicative of the flow. 
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□ SKEWED/SLANTED IMAGES 

□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE(S) or EXHIBIT(S) submitted ARE POOR QUALITY 

□ OTHER: . 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



BEST AVAILABLE IMAGES 




